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ABSTRACT 

The relationship between th2 applied wind stress and 
currents predicted by a primitive equation ccean circulation 
model was analyzed and comparei to theory and observations. 
Three cne-~year data sets wer2 examined using Fourier and 
rotary spectrum analysis techniques. The Fourier analysis 
revealed three spectral peaks in the predicted currents with 
none in the wind stress. Thes2 peaks correspond to synoptic 
Variability at low frequency, the inertial response at an 
intermediate frequency anda nonphysical response at high 
frequency, due to the finit2 time differencing procedure 
employed. This response at high frequency was «wo orders of 
magnitude smaller than the peaks at the synoptic and iner- 
tial periods. ie Cred’ emot2On was the same order of 
magnitude as the synoptic motion near the surface, but much 
weaker below. I+ was identified by the rotary spectrum, and 
it was Slightly shifted toward lower freqguencies in direct 
proportion to the time step used by the model. The Ekman 
moticn appeared to be restrict=d to the baroclinic response 
above 70 m. The time-varying "geostrophic" flow below 70 o 


was essentially barotropic. 
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I. INTRODUCTION 


Pecos tlGOnral.y, general ocean circulation models have 


been verified against the tamperature or density structure 
of the ocear. This approach was mainly =stablished because 
of the availability of ocean measurements. Temperature and 
density, through salinity values, are the ocean variables 
recorded over the largest area of the world's oceans today. 
As the number of bathythermograph recordings and ocean sta- 
sions increased, a pocl of temperature and salinity méasure2- 
nents became available, forming a general oceanic 
climatology which could be used to gauge 2n ocean ea = 
tion model's performance. Once reasonable confidence was 
established in the model's capability to reproduce known 
ocean thermal and/or density characteristics, sensitivity 
experiments cculd be conducted to test the effects of cer- 
tain changes in various model parameters on aifferent 
modeled features. 

Since the ocean is a vast sensible heat source or sink 
for the atmosphere, ocean circulation nodeling is not a 
peetomimee 2tselt, but it is also tisd to @mprove atmos- 
pheric modeling. This is a two way link because both fluids 


can be nodeled by the same general set of physical equations 
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and, depending on *he time scale examined, each forces the 
Owtler in thear overall jcint circulation. As a result, a 
better understanding cf the ocean's thermal structure, espe- 
cially its upper profile and effects which modify it, could 
lead to improvements in an ocsanically coupled atmospheric 
circulation model over the relatively data sparse areas of 
the worid's oceans. 

The objective of this thesis is to analyze the currents 
predicted by an ocean circulation model as part of a general 
model verification program. This analysis is expected to 
improve the ability to model and predict changes in the oce- 
an's motion, thermal fields and its response to atmospheric 
forcing on a weekly tc seasonal time frame. Besides their 
Own intrinsic importance, currents play an important roll in 
influencing «he thermal structure in the océan. In this 
respect, their contribution may be twofold: turbulent mixing 

nd large-scale advec“ion. Strong shear currents are res- 
Pomel LOr turoulent mixing <hrough convection, while the 
large-scale Ekman/geostrophic flow contributes +o advection. 
Besides improving its thermal profile, an ocean circula*ticn 
model which correctly depicts nass transport would have a 


greater cverall applicabilic than one which is just 
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concerned with displaying the thermal structure of a water 
column. 

One of the reasons currents have not been examined in 
verification studies earlier was the lack of available open 
eccean current data. But as part of NORPAX (North Pacific 
Experiment), @ data base of Sp2n ocean loig-*term currents is 
now becoming available, McNally (1981). AS Pal] of ches 
ongoing experiment, Gerald J McNally of Scripps Institute of 
Oceanography, La Jolla, California placed a number of drif<+- 


ers, both drceuged at various depths and undrouged, in th 


(D 


omen Pace tic Drift Current. Each buoy was equipped with a 
transmitter and was tracked for as long as six months by 
RAMS (Random Access Measurement System) onboard NIMBUS 6. 
Several results of this ongoing investigation are now avail- 
able. One, the apparent lack of difference between current 
observations at the surface and at 30 m, will serve asa 
qood test for examining «he current profile produced by an 
ocean circulation model. The observed statistical relation- 
Ship between winds and currents also nezds to be investi- 


gated from a theoreticai or modeling point of view. 
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The model evaluated was developed by Robert L Haney, 
Naval Postgraduate School, Monterey, California (Haney et 
al, 1978). It is a primitiv2 equation model applied to a 
baroclinic ocean in a rectangular basin. rt ‘Contas ns 20 
levels in the vertical with 8 in the upper 200 m. The nodel 
*s based on the hydrostatic and Boussin2sq approximations. 
A rigid-lid approximation is made at the surface and salin- 
*ty is neglected. The effects of subgrid or mésoscale 
moticn are parameterized by vertical and nacniinear horizon- 
zal eddy viscosities and conductivities. 

The model is driven by heating and stresses applied at 
the surface. The downward flux of heat and nomentum is 
determined by the following time-dependent variables: the 
predicted ocean surface temperature, the prescribed values 
of 6-hourly surface winds and monthly mean: air temperature, 
Specific shum-dity, fract:cmal clomid Seower awd ineoniag 
Shortwave radiation. These quantities are used *o calculate 
she shortwave and longwave radiation, sensible and laztent 


heat and surface wind str2ss. The buik aerodynamic 
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equations for determining wind stress from marine winds 
along with the primitive equations of motion, the continuity 
equation, the first law cf thermodymanics and the equation 
of state can be found in Haney, et al (1978). The governing 
solar radiation equations can also be found in this refer- 
ence with one exception. Anew approach as outlined in 
Haney (1980) allows scme of th2 solar radiation to panetrate 
to depths greater than 20 n rather than being totally 
absorbed in this layer as was the case with the old radia- 
tion equation. 

To handle the transport of heat and momentum which 
OCCCUIS 9h Scaies too small to be resolved by the grid spac- 
ing employed, parameterizations of these motions are made on 
scales which are discernible by the model. Since eddy flow 
is important, especially in was ern boundary currents and 
their extensions, this motion must not be neglected in an 
ocean circulation model. Until modsl grids are reduced to a 
size sufficiently small to allow eddy representation, vari- 
ous theories must be applied to approximate eddy motion. 
The Haney model uses the mixing length and momentum hypoth:- 
Sis of classical turbulence theory. This approach intro- 


duces coefficients Ome herrzo ntal eddy viscosity and 
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conductivity. Vertical diffusion is approximated through 
small coefficients of vertical eddy viscosity and conductiv- 
ity, see Haney, et al (1978). 

Active turbulent vertical mixing of heat and momentum in 
*he model is based on a parameterized "dymanic adjustment". 
The key to the adjustment is the local gradient Richardson 
number. This number is a non-dimensional measure of stabil- 
i=y as it relates the local hydrostatic stability or densix«y 
stratification to the square of the local shear. When “he 
water column is neutral, the local Richardson number takes 
on its critical value of 0.25. This values was theorstically 
established and is justified through observations presented 
in Thompson (1980). The adjustment consists of checking the 
water column at ali levels b2low the surface for dynamic 
Siisradpa Lc t Ye mers CONdit19n =Xists if Ehe local gradien= 
Richardson number computed between two adjacent model levels 
falls below its critical valué. Mien sa case, 6 Verru.e2. 


mixing of both temperature and velocity components between 


(Db 


the adjacent levels transpires in such a way that: (a) heat 
and momentum are conserved,  (b) the mixing ratios for heat 


and momentum are set equal ani (c) che adjusted gradient 


Richardsen number is set equal to its critical value. The 
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specifics of this adjustment procedure can be found in Ada- 
mec, et al (1981). It is through this parameterized mixing 
mechanism that «he model predicted currents affect the ver- 
tical turbulent mixing of both heat and momentun. This is 
of the shear of the mcdel predicted currents. 

Baroclinic (or shear) currents are handled through 
G2S@ct Calculations Erom the horizontal eguations of motion. 
The vertically averaged currents are defined by 4 strean- 
function derived by a vertical average of the horizontal 
Equations of motion and governed by the vorticity e2quation. 
These average currents are then added t9 the results of the 
primitive equations to yield the total currents. The model 
runs on a staggered grid with the temperature and stream- 
function gridpoints established at the model indices and 
momentum gridpoints offset by half 2 grid length. Time d:if- 
ferencing is accomplished by using ths leap-frog scheme 
with a Euler-backward scheme inserted everv 8 +imesteps +o 
reduce "time-splitting". 

The rectargular domain covered in tha model is in the 
Nova Nea cs £2C. Gowns SOS or eengr~cude in width, running 


poomretes 2? F to 1259 W and is 65° of latitude in length, 
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running from the equator to 659 N. This domain allows the 
model to simulate the circulation pattern of the interior 
Morth Pacific Ocean. The general surface flow includes «he 
equatorial flow pattern, the North Pacific gyre and the Sub- 
polar gyre systems. The complex flow around the Japanese 
and Aleutian Islands are not revealed because these intrud- 
ing land masses are not included in «he boundary conditions. 
The bottom is flat ata uniform depth of 4&4 kn. At the 
lateral boundaries, there is parfect insulation of both heat 
and momentum with zero velocity except at the equator where 
a free slip condition is impos2i. The bottom is also ther- 
mally insulated but a small stress is imposed through a sin- 


ple quadratic drag law, see Weatherly (1972). 
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Tit. MODEL OUTPUT 


A. MODEL INITIALIZATION AND SIMULATION 

The model simulation was primarily driven by 10 years 
(from 01 Jan 1969 through 31 Dec 1978) of 6-hourly marine 
wind analyses performed by FNOC (Fleet Numerical Oceanogra- 
phy Center, Monterey, California). These winds 2re the 
Mertabls Input to the bulk formula used to calculate the 
surface wind stress and the latent and sensible heat fluxes 
in the model. The Other e-omOspheric forcing variables: 
solar radiation, Giveud Gouonm, aic temperature and vapor 
pressure; are all prescribel by monthly climatological 
values. The sources of these values are given in Haney, ¢t 
ae (19750). 

The model started from an initial state with zero cur- 
rents and a temperature structure which decreased linearly 
with depth ranging from a valu equai <9 the local air temp- 
eracure at the surface to 29 C at the botton. It then ran 
for 3.5 years from 01 Jan 1974 through 390 Jun 1977 with the 
forcing being supplied by th2 FNOC wind analyses and the 


monthly atmospheric climatologies. HUrind = thas einscce! 
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phase, the model used 25 fixed levels in the vertical (12 in 
the upper 200 m) running on a coarse grii with a spacing of 
approximately 6° of longitude and 4° of latitude. 

After this initial "spin up" was complets, the model was 
run to simulate three decades of ocean circulation using the 
10 years of wind data available. The first 10 years of 
Simulation was run on a coarse grid (69x4°) and the succeed- 
ing 20 years were run on a fine grid (39x29). Bach =2ime a 
new decade of Simulation was Started, 2 smoothing was 
applied to the model variables (currents and temperatures) 
*=n order to reduce "model shock" when the wind forcing was 
shifted back to QO1 Jan 1969. This smoothing consisted of 
averaging the simulated variables over th= last 2 weeks of 
that particular decade and using these average values when 
the model was reinitialized to Jan 1969 to start the next 
decade 9f simula*ion. When the model was switched to the 
fine mesh grid, the grid siz> was cut in half and values 


were linearly interpolated tc the new gqridpoints. 


ie COARSE GRID OUTPUT 
The coarse grid used for the initial decade of ocean 
Simulation consisted of 17 points in both the longitudinal 


and latitudinel direction. AS a result, the grid spacing 
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WAS approximately 6° cf longitude and 4° 9£ lakitude. The 
msdeél was run with 20 levels in the vertical vice the 25 
levels used in the initialization phase. The levels used in 
this simulation (for both the coarse and fine grid) were as 
Be ome On og 31, 49, 71.5, 100, 136.5, 182.5, 240, 313, 
M@5-5, 222, 669, 855, 1090, 1385, 1760, 2235, 2830 and 3580 
m below the undisturbed surface. 

The variable coefficents of horizontal eddy viscosity 
2nd conductivity as displayed in Haney, et al (1978) come 
from the following equation - 

A= A.f{i+xXi7s| (as) | (i) 

where $ is the relative vortisity of th= vertically aver- 
aged current and AS represents a constant which is a meéa- 
sure of the grid size. In the present model ~~ = 4 x 10 
sec/cm. For the coarse grid, the baskground value for 
momentum is A, = 5 x 10° cm /s2c and for heat is 1x 10° cn 
JSEC. The vertical eddy viscosity is %, = 5 cm /sec and 
SOmauer  Viry 25 K, = 0.5 cm /s32. The coarse grii is run 
using a 6 hour time stéep. 

During this first decade, one sampl2 of model output for 
Siiwmelevers =n the vertical 2 a single gridpoint was 


extracted for analysis. mis JELAdpSint shoosen was located 
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NORPAX SATELLITE TRACKED 
DRIFT EUOY TRAJECTORIES 
| Segt. 10,1976 ~ Aug, 21,1977 





at 42.669 N and 150.319 W. Ses Fig. 1. The analysis period 
went from O17 Jul 1976 00Z through 30 Jun 1977 182. Ths 
time frame and location was choosen for a preliminary analy- 
sis cf the model's reponse because it corresponds tc the 
space and time domain of the drifter data collected by G. J. 
McNally for NORPAX and is in the major strom track across 
Pe WemN@iaci PaCa £1C. Tt is possible that the initialization 


of she nodel and the transition to the first decade affected 
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“he quality cf the output used in the analysis. But since 
the model had run through 8.5 years of continuous simulation 
before the output was extractsd, and because *he analysis 
was Only conducted decwn to 313 m, it was felt that any 
influence of the initial conditions would have long since 


been dispersed from the current structure. 


on FINE GRID OUTPUT 

The fine grid representation was est2blished to improve 
the horizontal resolution of the predicted variables and the 
*nertial respense in the model. This was accomplished by 
reducing both the grid spacing and the tine step used in the 
finite differencing scheme. 

The fine mesh version us2s 33 gridpoints in both the 
longitudinal and latitudinal directions. This reduces the 
grid spacing to approximateiy 39° of longitude and 2° cf 
latitude or twice the resclution of the coarse grid. Once 
again, a staggered grid is us2d with the temperature values 
determined at the model indices and the momentum values off- 
Sec bP halt ¢ grid léngth. The time st=p used in integra- 
tion was also halved «to 3 hours t> maintain linear 


computational stabiliy. 
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The only cther change is in th2 magnitude of. the hori- 
zontal eddy viscosity and conductivity cosfficients. These 
constants are reduced by a power of ten t9 a value of A, = 5 
x 10° cm” /sec for viscosity and 1 x 10 cn’ /sec for conduc- 
tivity. The governing equation may be found in the previous 
section covering the coarse grid. 

All the other model parameters and criteria are the same 
as the coarse grid version. The same 20 levels are used for 
vertical resolution. The relatively small coefficients of 
eddy viscosity and conductivity in the vertical have also 
been left unaltered. 

The fine mesh Simulation of ocean circulation consis*ed 
cf «wo decades of model execution. After the two week aver- 
age from the end of December 1978 was applied *o *hz2 coarse 
grid at the conclusion of the first decad2 modeled, the grid 
length was converted to the fine mesh Length. A linear 
interpolation of the averaged u and v currents and tempera- 
tures was performed to assign values to the newly created 
intervening gridpoints. Mine second décaede of current sinu- 
lation was chen conducted reusing the 10 years of marine 
wind analyses supplied by FNOC. T> wus2 2 33-hour *ime step 


=f te Tene mesh run, a linear interpolation of the 6-hourly 
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synoptic winds was performed. The other atmospheric forcing 
parameters still came from th? same monthly climatology as 
used in the coarse grid run. 

When the second decade (first with the finer resolution) 
was completed, a third decade was run. Ince again, the two 
week average smocthing procedure was employed en the ocean 
variables between decades before the winis were reinitial- 
ized to 01 Jan 1969. It was from this third decade of simu- 
lation that the two fine mesh samples of model output for 
data analysis. 

The firs+ fine grid data set ran from 01 Jul 1969 032 
tEncowoh OT Jul 1970°00Z. It consisted of a single gridpoint 
with all the values in the vertical beiag extracted as in 
=he coarse grid case. The grid point extracted in this case 
was located at 43.679 N and 148.919 W. Se cr ques The 
analysis of this data allowed an interannual comparison in 
studying the model's fine mesh performanc:. 

The final data set matched *he timeframe and general 
location of the coarse grid set. It was located at 41.64° N 
arammiont. 72° WW. See Fig. 1. This set was studied in «he 
areatest detail in an attempt <> understand the model's tes- 
ponse and compare i* to some of the preliminary results 


available from the NOEPAX drifter data. 
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A. ANALYSIS PROCEDURES 

Each model data set extracted from the simulation was 
analyzed using a new Biomedical Computer Program, P-series 
Subroutine, BDMP1T - Univariat2 and Bivariate Spectral Ana- 
lysis. This program package became available in th begin- 
Reng of 17962. It uses cosine transforms ina fourier 
analysis approach to reduce input data into its spectral 
components at all the analyzable frequensies, from 0.0 to 
0.5 (the Nyquist frequency). The subroutine's output fre- 
quencies are dimensicned per unit data time step and asa 
result must be converted +o be meaningful to the particular 
data set analyzed. The spectral packag2 runs its analysis 
using three different bandwidths: narrow, medium and wide. 
These bandwidths may either b= individually specified or 
assigned by the program if left to the default mode. In the 
default mode, the bandwidth size is dependent on the number 
of values in the input data set. Since this study was only 
concerned with the large-scale trends displayed by the Haney 


general circulation model, only the widest bandwidth results 


Z> 





were utilized. These results ware also the smoothest curves 
produced by the subroutine and are statistically the mest 
stable. For a complete description of this program package, 
refer to Engeiman, etal (1981). This spectral option was 
chosen over the other program libraries available at the 
Naval Postgraduate School because of the ease in interpre<- 
ing its output, variety of analysis »yptions available, 
detailed description in *he user's manu2l and clarity in 
diagnostic messages produced. 

A detailed analysis of the fine mesh 1976-77 data set 
was performed. It involved a description of the temperature 
and stability profile of the simulated upper océan over the 
annual cycle, a spectrum analysis of the predicted currents, 
a lowpass spectral analysis to 2xamine th: model's wind/cur- 
rent relaticn and a rotary sp2ctrum analysis *o examine its 
ineztial respcnse. The lowpass analysis was a two stp pre- 
cess. Paar St), *he model curren= data was run *hrough an 
eleven-point lowpass filter +¢9 remove the energy above the 
inertial frequency -as determined by *he original specctrun 
analysis. The weights used by this lowpass filter were: 
emcee o2, 0.0799, 0.1231, 0.1561, 0.1684, 0.1561, 


Oviizesen Ow0799, 0.0402, 0.0166. After lowpass filtering, 
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the 3-hourly data was resampled twice a day and spectrally 
analyzed once again. This procedure resulted in increasin 

the discernibie analysis range produced by the Biomed spect- 
rum analysis package. With a sampling interval of 12 hours, 
the analysis period was increased to 1 t5 20 days as con- 
pared to the 6 to 120 hour range depicted by the 3-hour time 
step data of the general study. The relationship between 
the input wind stress and near surface currents was then 
examined at this extended frequency range. 

The rotary spectrum analysis produced 2 decomposition of 
the model spectral currents ints their cyclonic and anticy- 
clonic components. This procedure utilized the autospectrum 
of the u and vcurrent components and their quadrature 
spectrum at each of the model levels studied. The result of 
Saeecmemsyscas iS a ratio of the cyclonic to danzicyclonic 
circulation which can then be compared ‘> the theoretical 
ratio using the analyzed frequencies and *he inertial fre- 
quency. Mecerd ing» to theery (Mecers, 1973), the value of 
<his ratio approaches one at the very highest and lowest 
frequencies (compared to the inertial frequency) and is zero 
at the inertial frequency. Because of this, the modei'ts 


inertial frequency wili be indicated by 2 minimum in the 


H 


atic. The equation used is - 
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R (s) = S(s) (2 ) 


GF} ~ 


the cyclonic spectrum: S(s) 0.5 (Guu (s) +Gvv (s)+2Quv (s) ) 


anticyclonic spectrum: S(-s) = 0.5 (Guu(s) +Gvv (s) -2Quv(s) ) 


where s20 is a given frequency, f is the inertial frequency, 
Guu(s) and Gvv(s) are the autospectra of u andv respec- 
tively and Quvfs) is the imayinary part (quadrature spect- 
rum) of the cross spectrum between u and v. 

The results of these analysis procedures will be 
described and interpreted in the following sections. The 
first data set treated is the fine mesh 1976-77 data. A 


generai description cf the relevant model characteristics 


wili be presented, followei by the spectrum analysis 
e=sulis. The model's inertial frequency will then be dis- 
cussed, eoeowed DY SUpperting results from the rotary 
spectrun. Then the wind/curr2nt relation will be presented 


using the lowpass filtered analysis. 
Once the analysis of the 1976-77 fins mesh data is pre- 
sented, a comparison will be made with ths 1969-70 fine mesh 


data. Aiter this comparison, the 1976~77 data set will b=: 
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compared to its coarse grid counterpart. The results sug- 
gest there is little difference between the two different 
years of fine mesh data but some important differences 


between the coarse and fine grid data sets exist. 


Bae CONE MESH RESULTS 
1. General Description of Some Output Parameters 

Each of the data sets was analyzed at four model 
levels: level 1 (5 m, called the surface level in the res- 
luts), level 3 (31 m), level 5 (71.5 m) and level 10 (313 
Br). These particular levels were chosen to discern the 
model's characteristics through its approximate mixed layer 
representation, in the transition zone of thermocline below 
*he surface mixed layer and in ‘the interior of the medeéel. 
This dscompesition of the watar column at the gridpoin:t 
selected is supported by the temperature, stability and cur- 
rent structure produced by the model. 

mee = cemperatire structure during 1976-77 is dis- 
played in Fig. 2. It showS a typical mixed layer structure 
going from a very warm, shallow (approxinately 10 nm) mixed 
layer in the summer months of June, July and August 720 4 
deep isothermal layer in th winter and early spring, 


approxinately February through April. The mixed layer 
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Figure 2. re Mesh Temperature Structurs (JUL 76 - JUN 


Fi 
77 
extends te at least 30 m from October through May supporting 
oe selection of 30 mas a good depth tor approximating he 
Central part of the mixed layer. The selection of depth 
m5 eas themkean transition zoné is supported by the loca- 
tion of seasonal thermocline, which is below 71.5 m from 
August through January. The seasonal tharmocline does not 
reach <hat depth in the summer and vanishes in the winter, 
as expected. At 313 m, the tamperature is relatively con- 
Slan= sirowenocut the year, placing it below the surface 
thermal effects and indicating that it represents “he common 


water in the interior of the model ocean. 
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1976 1977 
Figure 3. Richardson Number based on Total Currents for 
ite Mesh Model (JUL 76 - JUN 77). 

The Gymanic stability of the ws#ater column, as 
tTepresented by *he Richardson number, also supports the 
thermal structure breakdown criterion selected. Egy. 3 
shows the gradient Richardson number (Ri) jin increments of 


0.5 based on the «otal current, where - 
(3 
Ri <332/[(3 33) + (33 in “) 


and X is the thermal expansion coefficient, 9g is the gravi- 
tetional acceleration and Tis the temperature. The contour 


Dieemuats:  Pplorcted represent &2=0.5. The smallest value 
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contour of Ri=0.5 indicates that the turbulent portion of 
the water column is abcve this line. The critical value of 
Ri=0.25 (where the water column becomes neutral) lies in 
*his upper region and marks the base of this surface turbu- 
Noet Structure. As the Richardson number increases, the 
stability of the water column increases. The largest value 
contoured represents Ri=5.0. Tsing the iepth of the criti- 
cal Richardson number as an indication of the <“ransition 
zone below the well mixed layer, it is noted thet chs 
level is deeper than 71.5 m in the winter and shallower in 
summer. This shows the selection of 71.5 m as an average 
for the transition zone was reasonable. The importance of 
ene Shear of the inertial motson in determining «he stabil- 
ty can be displayed using a plot (Fig. 4) of the gradient 
Fichardson number calculated from averaged currents with the 
interial frequency removec. As can be seen, the unstable 
zone does not penetrate as deeply as it does in Fidq. 3. The 
G@S—peh Of turbulent maxing that would exis: in the absence cf 
inertial motion is somewhat less accordirg =o this compari- 
Seu Thus, these motions make an important, though perhaps 
hot Gummadi, COntE. Dt: Oon to zche turbulent méexing in this 


model. 
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The above result is also supported by the current 
profile produced by the model. As seen in both the east- 
west component, Fig. 5, and the north-south component, get ete 
6, the major current fluctuations occur above approximately 
30 a. Below this depth, the currents show little change in 
the vertical ard are reduced in magnitude. Fig. 6 displays 
a periodic dizection reversal bslow about 70 m, as demons- 
rated by the heavy vertical line representing a zero 
north-south velocity. Diese pat csrh depicts a bvasteally 
barotropic response below the mixed layer which is also seen 


in the spectrum analysis presented in the next section. 


Northward Current (cm/sec) 
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A spectral analysis of both the input wind stress 
and the currents produced by the model was made. The band- 
width used was 0.2816 cycles/day (cpd). The wind stress @s 
seen in Fig. 7 was analyzed to determine if a spectral peak 
was present which could be transmitted to the ocean and show 
up in the current Simulation. With no peak found, it can be 
safely assumed that any spectral peak in the currents is not 
a direct result of the wind stress forcing and must be due 
zo some oceanic "natural selection" process. Figs. 8 and 9 
display the power spectrum cf the u and v current compo- 
nents, respectively. Three spectral peaks can be seen. The 
ficst is at relatively low frequency and represents the 
transient flow oroduced by atmospheric synoptic forcing in 
“he model. This peak will b2 examined in the wind/curren= 
relation section. The beweetongocmorcton in the interior of 
=he modeled water coiumn displayed in Fig. 6 is represented 
DY £RES portion of the speéctrun. At frzquéncies less than 
0.30 cpd, the energy levels a+ level 5 (71.5 m) and levei 10 
(313 m) are the same indicating that most of she energy in 
<his frequency range is confined above 71.5 m and the model 


response is basically barotropic below this depth. The 
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upper portion of the ocean displays baroclinic conditions in 
contrast to this barotropic response in the model interior. 
At the lowest frequencies (<0.1cpd), the energy spectra in u 
and v are of the order of 2.653 x 10 and 1.761 x 10 at the 
surface and 5.679 x 10 and 3.799 x 10 cm* /sec” at 71.5 o 
and kelow, respectively. This indicates RMS currents at 
these frequencies of about 27.3 and 22.3 cm/sec at the sur- 
face, with 4.0 and 3.3 cm/sec at 71.5 m and below. 

The dominant peak in the spectrum is centered in the 
e570 =) 1.0256 cpd range. Tts amplitude is of the same 
order of magritude as the low frequency peak in the upper 
three levels analyzed. This peak is also evident at 313 a, 
but here its amplitude is much smaller than that of the iow 
frequency motion. Since the inartial frequency represents a 


high energy level in the ocean compared to the fréquencies 


around it, -+ is concluded that this peak is the nodel's 
inertial response. This response is discussed in greater 
detail in the next section. The last peak visible is very 


small in amplitude, lecated at relatively high frequency and 
appears only in the urper levels studied. This peak is out- 
side thea objectives of this study and is only investigated 


because of the marked effect ln the rotary spectrum 


36 





associated with it. It appears at a frequency corresponding 
to a period roughly 2.5 times the model time step. It is 
PorenwcnasSenagn treaquency p@ak is an artifact of the finite 
difference scheme and not the result of any real physical 
phenomenon. Because it is at least two orders of magnitude 
Smaller than the other significanz features in the model, 
the overall effect of this peak on the nodel's performance 


should be small. 


Using the corlolis parameter calculated from the 
ieee tude Gf the gridpcint as input, the inertial period is 
laweez69 6 pd. The model's inertial response lies in the 
02.9576 - 1.0256 cpd range according to the spectral results 
and is further refined by the rotary spectrum to 1.0256 cpd. 
This shift toward lower frequencies appears to be du2 to the 
<rapezoidal implicit time diffarencing scheme employed by 
the model and is directly dependent on <=he length of the 
time step used (see comparison with coarse grid model 
below). 

AS stated in the previous section, the inertial 
moticn in the ocean will be associated with a dominant spec- 


tral peak when the curren*s ar2 studied over a year. This 
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uS Gilewcre Bene sendency £LOr intermittent ocean motion at the 
fnertial frequency. The signature of inertial motion in the 
model can be seen in the spectrum analysis through 313 nm. 
The energy associated with this motion appears somewhat 
masked in the surface level, is dominant in the mid levels 
and falls off at 313 m The spectral energy associated with 
half the height of each inertial peak at 1.0256 cpd is 
listed in the table below along with the frequency band over 


which this energy is distributed. 


TABL2 I 

Ne en TARE OSA ae at an a ee 

U-COMPONENT V-COMPONENT 

LeveL ENERGY BANDWID?H ENERGY BANDWIDTH 
sfc (5m) eos x 10° #48©60n238 3.981 x 10 0.446 
L93 (31m) 1.047 x 10 0.349 7.943 x 10° 0.365 
Poot? ieom) 7.096 x 10 0.360 7.244 0.376 
L10(313m) 7.244 x 10° 0.287 6.607 x 10 0.303 


The Spectral phase information produced when the u and v 
components were analyzed shows that the inertial motion in 
=he model wae anticyclonic. The v-coOmponent leads the 
u-component by approximately 9099 (from 89.79 at the surface 
+9 90.69 at 313 m) For all four levels studied, as expected 


zteom theory. 
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The rotary spectrum confirms this spectral peak as 
the model's inertial responses. Gre eeneoey sbehand this ana- 
lysis approach (Mooers, 1973) dictates that the ratio 
between the cyclonic and anticyclonic flow components ina 
water column go to zero at tha inertial frequency where the 
motion is purely anticyclonic. The results of the rotary 
spectrum calculated from equation (2) f sr each level ana- 
iyzed are depicted in Figs. 10 - 13 respectively. The 
theoretical curve calculated from equation (3) Weng she 
actual inertial frequency calculated from the gridpoint's 
coriolis parameter is also plotted for comparison. ic Can 
be cencluded from these curves that the nodel's best perfor- 
Bemee, e@ecord’ng to the supporting <heary for “the rotary 
spectrum, is in the frequency range of 2pproximately 0.3 - 
Dace pa. From the tables of analyzed frequencies used +o 
produce these curves, it can be seen that the minimum in «he 
Gese@uicted FTacio iS at 7.0256 ¢pod which corresponds to the 
inertial frequency determined from the spectral analysis 
covered in “he previcus section. The tables used in the 
rotary spectrum may be found in Appendix A. Since most of 
the model's energy at each level is contained at frequencies 


lower Shan 2.2 cpd, its response is generally good when 
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examined by this methcd which is based 91 inertial-internal 


wave theory, £<O<N. 
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Figure 10. Fine Mesh Rotary Spectrum for Surface (5 mn), JUL 
76 - JUN 77. 


According to the rotary spectrum, the high frequency 
peak visible in the spectrum analysis (near 3.2 cpd) appears 
+>) be asscciated with strong cyclonic motion at all levels 
Studied. But the RMS velocities at this frequency are: 1.8, 
Seon eame. > cmysec at the suetace, 31 m and 71.5 nn, tespec- 
Sey There is no peak evident at 313 on. As alluded «+o 


earlier, this motion, having a2 period near 24t, may be the 
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esh Rotory Spectrum for Level 03 (31 n), 
JUN 77. 


eM 
"oe = 


pagire ty, Fin 
JUL 
result of time splitting in the leapfrog time differencing 
scheme. The model's performance can b2 improved if «his 
peak is removed; but because of the very small 2mount of 
energy associated with it, its impact on the modei'ts general 
ccean circulation charactéristics should be minimal. 
4. Wind/current Relation 

The aim cf this part of «he study is to inspect the 
modél's performance in regards to Ekman theory (wind driven 
flow to the right of the surface stress in the Northern 


Hemisphere). Accordingiy, the osutput current structure will 
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Figure 12. Fine Mesh Rotary Spectrum for Level 05 (71.5 m), 
JUL 76 - JUN 77. 


be examined where «there is high coherenc® between the cur- 
rent and its orthogonal wind stress. The point at which 
this relation breaks down in going towari high freyuencies 
=ndicates the limit of steady Ekman response and the transi- 
tion in governing physics to time idependen~ Ekman response 
=o the applied wind stress. Phas transition wiil be sup- 
ported by an increase in the coherence with the parallel 
wind stress. It is also expected that the coherence between 
thee GMrrcent end its orthegonal wind stress will rapidly 


decrease with depth below the model's mixed layer, 2ssuming 
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Pugire 13. Bee Fes) RO Laay Spectrum ror Level 10 (313 n), 
this level is a good first order approxination of the Ekman 
@epth. Meee aepth, the Er2cticnal effect of «he applied 
wind stress falls off to4 percent of its surfaces value. 
Below this depth, the geostrophic componen= of *he current 
should be removed from the flow before an evaluation of any 
relation between the applied wind stress and currents is 
made. 

The results of the coherence study between the cur- 


Tents at the selected model depths and «he applies wind 


Smecss using ~he lowpass filtered date running from 90.05 - 
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1.0 ecpd generally ccnformed to the hypothesis presented 
above. The coherence was generally high between the u-con- 
ponent of the current and *he y-component of the wind 
stress. This coherence pattern was mimicked by the v-cur- 
rent and x-wind stress except at lower coherence values. 
As expected, the coherence falls off with depth with a 
marked decrease between level 3 (31 m) and level 5 (71.5 m). 
Recall that level 3 provided a satisfactory indication of 
the base of the mixed layer when averaged over the entire 
year of data and level 5 provided a satisfactory indication 
of the transition zone or seasonal thermecline below the 
mixed layer. According to th? coherence results, it main- 
*ains a coherence of less than 0.4 over most of the fre- 
quency range examined and ther2fore appears to be primarily 
below the Ekman depth in «he model. Because cf the strong 
esascnal signal, a spectrum analysis of seasonally strati- 
fied data may be beneficial in the future. 

The coherence values between the surface u-curren= 


and y-rwind stress ranged from 2 high of 0.966 to a low of 


Ozes3 The coherence is over 0.90 in the frequency range 
P=Eemmenoe—= O.2625 Cpd and finally drops below 0.80 at 0.387 
cepa ae level WP (5 °m). Meco tenes DOInRt, +he coherence 
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continues to fall through the rest of the analyzed frequency 
range. This overall pattern is reflected in the level 3 (31 
m) analysis, only with lower coherence values. The maxinun 
value is reduced to 0.853 with a minimum of 0.090. The 
values stay above 0.75 from 0.0 - 0.516 cpd and finally drop 
below 0.65 at 0.90 cpd. This trend completely disappears at 
level 5 (71.5 m). The coheresnc2s are all below 0.375 except 
at two peaks - one centered at 0.57 cpd, the other at 0.856 
God. These peaks increase to coherences of 0.650 and 0.575 
respectively, and may reflect Ekman pumping/suction. They 
ace also visible at level 10 (313 nm). From this analysis, 
it appears that the model currents follow Ekman theory at 
the surface and level 3 at frequencies as high as 0.90 cpd. 
At higher frequencies and at or below 71a, there is little 
meaningful reliationship between the steady Ekman response 
and the currents. 

The phase information produced by the ccerrelation 
study between the wind stress and tha nodel currents was 
examined to investigate the deflection angle produced by the 
model's Ekman response. Only the phase angles at F = 0.0 
cycles/day corresponding *o the annual signal were used, but 


the results are valid inthe mod2l for time scales nuch 
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greater than a day. Graphical solution techniques wers 
employed using the spectral energy at 0.9 cpd as *he magni- 
tude of each component to determine the deflection angle 
between the total wind stress and the currents at the sur- 
face through level 5 (71.5 mn). There is a lot of variation 
in the results because, though the coherence is high between 
the currents and their orthogonal wind stress (approximately 
0.88 above 31 m), it is low in the cross-correlation between 
both the current components and th wind stress components 
themselves (abou~ 0.18). 

In the three levels examined, the phase was 0.09 
between the u and v-currents; between u and the y-wind 
stress, and between the x-str2ss andthe y-stress. The 
faaeccummens X-Wind stress were 180° su* of phase. ie 
v-current leads both the u-current and x-wind stress. The 
y7>wind stress leads the u-current. When the wind and cur- 
rent components were plotted, some interesting results were 
noted which confirmed the conclusions reached above. The 
surface current was deflected roughly 83.09 to the right of 
the wind stress, with the current at 31 m deflected 14.59 
Pita ene= for a tonal of 97559 to the raght of the wind. ole 


the current at 71.5 m was approximately 24.4° to the left cf 





tae SOM Ceeegnt or 73.19 t> the right.»sf the wind. The 
expected Ekman spiral representation breaks down at 71.5 nm. 
Bux, recall the coherence study revealed very littl? corre- 
laticn between the wind stress and the currents at this 
depth. This may be a further indication that the Ekman 
depth in the mod¢l is above 71.5 mo. It is felt this depar- 
ture from Ekman theory is due to the inclusion of the geos- 
trophic currents in this study. The geostrophic currents 
should be remeved before any concrete conclusion is reached, 
but it appears from this analysis that the open ocean cur- 
rents are consistent with Ekman theory at low frequency. 
The deflection angles near 909 are reasonable because the 
model assumes a well mixed upper layer (1.3. a uniform ,vel- 
Gerty, Tet unzfoma eddy viscosity) and therefore the cur- 
rents €quate to the mass *ranspor= at each level. They do 
rot indicate a surface current which would be deflected 45° 


mo thesriight of the wind according t> Ekman <=heory. 


Cr FINE MESH INTERANNUAL COMPARISON 

The major finding in the interannual comparison was that 
Ptemwrsp2=evral tceatures establishei in the "1976-77 data 
remained virtually unchanged in the 1969-70 data. Just the 
magnitude of the spectral intensities were reduced in 


NO 9=70. 
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Once again, the input wind stress was analyzed to deter- 
mine if a spectral peak was prasent which could be reflected 
in the model currents, see Fig. 14. When no peak was found, 
@ comparison was made with the 1976-77 wind data. The pat- 
tern was similar, but the spectral energy associated with 
the 1969-70 winds was consistently lower than its 1976-77 
counterpart (except at frequencies below 0.2 cpd where the 
x-cOmponents were equal). This demonstrated an intensity 
difference in the low frequency winds from one year <o 
another. 

The reduced intensity in the spectral energy associated 
with the winds is reflected in the predicted model currents. 
The three spectral peaks associated with Ekman flow, the 
inertial frequency and the time differencing error at high 
Peeaiency appear in bcth the @ Curten*, Fig. 15, and vecug— 
rent, Fig. 16, analyses. The inertial frequency is still 
shifted toward lower values compared to the frequency calcu- 
lated from the ccriolis parameter associated with the grid- 
pount's Latzatude. The calculated inertial frequency is 
1.3806 cpd while the model's inertial response is at 1.0256 
cod, as it was in the 1976-77 rasults. The spectral energy 


is reduced atall levels and is generally confined in a 
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narrower frequency band than its 1976-77 scounterpart. This 
result is evident when the energy from half the height cf 
the spectral peak at 1.0256 cpd and the bandwidth over which 
-¢+ is distributed displayed in Table II is compared to those 
in Table I. 

TABLE If 


Half Power Spectral Energy (cm* /sec* ) and Bandwidth ios) 
for Model Inertial Mction for Fine Mesh (JUL 69 - JU 0). 


U-COMFONENT V-COMPONENT 
LEVEL ENERGY BANDWIDIH ENERGY BANDWIDTH 
sfc (5m) 2.754 x 10 0.195 1.995 x 10 0.338 
L03 (31m) 6.607 x 10 | OMe TH 5.754 x 10 0.337 
LO5 (71.5m) 3.020 0.410 eet 0.411 
maonisiem) 30631 2% 10) OL246 3.631 x 10° 0.223 


The high frequency peak is at the same frequency as the 
peak obtained in the 1976-77 data set study. This face is 
consistent with the hypothesis that it is due to a time dif- 
Perencing Error. The low frequency peak corresponding «=o 
Ekman flow is virtually unchanged in the two data set ana- 
lyses. The RMS u and v currents are respectively: 23.5 and 
20.6 cm/sec at the surface and 4.0 and 3.8 cm/sec at 71.5 om 
ard below. 

These Pi: indicate thac though the energy associated 
with the model depicted currents may change from year to 
year, its response +o the forsing paramters appears to be 


consistent, at least in the two years analyzed. 
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D. COARSE GRID COMPARISON 

When the 1976-77 3-hourly fine mesh data were compared 
so its 6-hourly coarse grid counterpart, some interesting 
conditions were observed. The 2xpected spectral pattern was 
repeated in the coarse grid analysis; but the siaqnal associ- 
ated with the spectral peaks was stornger compared to «he 
background spectral energy than it was in the fine mesh 
spectral results. 

The spectral peak associated with the modeled inertial 
respense was the most dominant feature visible in both the u 
and v current component analyses, Figs. 17 and 18, trespec- 
tively. This model characteristic is ‘the same as the fine 
mesh response; but it is shiftad away from the true inertial 
vericd toward even lcwer frequencies than the shift in <he 
fine mesh. Recall that the calculated inertial frequency 
for the fine mesh gridpoint was 1.3289 cpd and the model's 
response was at 1.0256 cpd. In the coars? grid representa- 
tion, the caiculated inertial frequency is 1.3549 oecpd and 
*h2 model's response is at 0.7220 cpd. The results of ches 
analysis suggest that a timestep »9f one (1) hour may be 
required to simulate the inertial motion a* much more nearly 


she proper frequency. The energy associated with half the 





height of the spectral peak at 0.722 cpd and the bandwidth 


over which it is distributed is displayed in Table IIT. 


TABLE III 
Half Power Spectral Fnerg (cn /sec’) and Bandwidth (cpd 
for Model Inertial Motion For coarse Grid (JUL 76 - 05 oh 
U-COMPONENT V-COMPONENT 
LEVEL ENERGY BANDWIDTH ENERGY BANDWIDTH 
sfc (5m) 7.079 x 10 0.100 5.888 x 10 0.123 
LO 3 (31m) 8.511 x 10 0.156 7s 10 nO oe 
L05 (71.5m) 8.710 0.168 8.710 0.123 
L10(313m) 5.888 x 10. 0.120 6.507 x 10. 0.108 


When the above values are compared to their fine mesh coun- 
terparts in Table I, it is noted that the anergy is slightly 
lower but is ccnfined in a narrower band in the coarse grid 
model. 

The high frequency peak found at roughly 2.5 times the 
time step in the fine mesh run Was at nearly 3 times the 


+ime step in the coarse grid run. Plus, it was amplified in 


}-4 


the interior layers and became a dominant feature at 71.5 ao. 
Because it was located at roughly the sam? point relative +o 
the time step and its maqnitudis was dependent oon the model 
resolution (including the timestep), this high frequency 
peak seems to b= due to the ‘time splitting characteristics 


of the leapfrog time differencing scheme used in the model. 


One Yappecach to test this hypothesis would be to run “he 
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model using the Euler-backward-step every iteration rather 
than every eight steps, to ramove the "time-splitting" in 
the model soultion. If this hypothesis is cerrect, this 
methcd would remove this peak from the spectral profile. 

The low frequency motion is ccnsistent with the motion 
display2d in *he fine mesh simulation. Th2 upper portion of 
the model ocean shows a baroclinic respo1is¢ with barotropic 
motion below 71.5 m. The u.andv power spectrum respec- 
tively display RMS currents of the order of 23.1 and 19.4 
cm/sec, respectively, at the surface and 2.5 and 1.2 cm/sec, 
respectively, at 71.5 m and b2low. These values are some- 
what lower shan those obtained from the fine mesh solution 
which is consistent with the reduced amplitudes depicted in 
the coarse grid inertial response. 

The result of doubling the resolution is multifaceted. 
wre trequency OL the inertial motion is closer to «he calcu- 
lated value in the fine mesh model, but it is smeared over a 
broader frequency range. Pie =sserec. Of the hagh frequency 
Wreme=-splitcting™ is reduced by increasing the resolution. 
The Ekman respcnse is Similar in both cases but the energy 


levels are sligntly higher in zhe fine mesh simulation. 
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The Haney ocean general circulation model consistently 
produces three spectral peaks when analyzed over a one-year 
period. These peaks are located at the lowest discernible 
frequencies, at roughly a mid-range vali2 and atthe high 
frequency end. They respectively correspond to the model's 
Ekman r2sponseé to atmospheric synoptic forcing, inertial 
motion and a possible numerical error due to the finite dif- 
ferencing method used. This overall pattern retained its 
structure when an interannual comparison was conducted 
uSing the fine mesh Simulation. Qnly the magnitude of the 
spectral energy changed in both the analyzed input wind 
stress and the predicted currants. This demonstrates «he 
variation in both atmospheric forcing and model response 
from year-to-year. When a@ comparison between the coarse 
grid and the fine grid was mad2 over the same timeframe, the 
three-~peak spectral pattern was repeated but amplified, 
shifted toward lower frequencies and concentrated over a 


hazrrower frequency band. 
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The inertial moticn in the model Simulation was shifted 
toward lower frequencies in both the fine and coarse grid 
representations. The frequency of the model's inertial 
motion was verified by a rotary spectrum analysis as deter- 
mined by the minimum value of the ratio sf the cyclonic to 
anticyclonic energy. In the 1976-77 fine grid simulation, 
the inertial motion was shiftei approximately six hours fron 
a true period of 18.06 h to 23.40 h in the model. In the 
i969=705daca the shift was from 17.38 h to 23.40 h. Rheic 
shift was even greater in the soarse grid resolution, but 
the inertial energy was confined to anarrower frequency 
range. The coarse grid response was shifted about 15.5 h 
Sa0m taecaue period of 17.71 h to 33.24 h. [has oShlte..2s 
Dartially due to the stability criteria of ‘the finite dif- 
ferencing scheme, which depends on fAt. With a 3-hour time 
step in the fine mesh representation and an inertial period 
S=weppssxiferely 18 hours, this accounts for roughly 3/18 or 
a 162m shift. A 6-hour time step in the coarse grid 
aceoumes for abomt 6/18 or a 33% error. 

The modeled inertial motion do¢s display some character- 
istics observed in «he ocean. The spectral energy associ- 


ated with this motion is high in the surface leveis and 
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falls off with depth. One governing condition appears to be 
whether the analyzed level is above or below the seasonal 
thermocline. The inertial enargy in the currents falls off 
rapidly below the surface layer. Such a decrease in iner- 
tial motion energy with depth was observed off Vancouver 
Island, British Columbia (Thomson and Huggett, 1981). 
Through the Richardson number study depicted in Figs. 3 and 
4, it was demonstrated that, though the contribution may not 
be large, the inertial motion does add to the turbulent mix- 
ing in ‘the model. The depth of the turbulent layer is 
reduced if the inertial motion was filtared out of the total 
model predicted currents when calculating the gradient Rich- 
ardson number. 

The low frequency responsé is consistent with Ekman 
AGW This was confirmed through a correlation study using 
lowpass fiitered data. This study showed 2 high correlation 
(a coherence on the crder of 0.88) between the predicted 
currents and the orthcgonal wind stress component down to 31 
M. The next level analyzed, 71.5 m, showed considerably 
less correlation (a ccherence of 0.37) with the orthogonal 
winds. This matches the preliminary findings from the 


NORPAX drifter data which demonstrate little change in 
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current observations over the first 30 m of the analyzed 
water column. The phase information from the correlation 
study indicated a mass transport to the right of the wind 
stress as expected freom Ekman Theory. The deflection angls 
was roughly 83° at the surface increasing to 97.59 at 31 @. 
But it decreased to 73.1° at 71.5 ma, indicatirg some other 
physics (possibly geostrophic flow) was also involved. 

To further delineate the model's current response to 
wind stress, the predicted currents could be divided into 
s2asonal groups. This would allow a better study of the 
Currents in relation to the nixed layer than can be done 
over an entire annual cycle. If all the levels were ana- 
lyzed vice a few selected ones, the bottom of the surface 
response could be located and compared to the NORPAX drifter 
data and other current measurenents in the ocear. Since «he 
geostrophic current components can be calculated exactly by 
ome emedel, <=hey could be subtracted from the total surrents 
leaving the ageostrephic current. This ageostrophic field 
could then be analyzed to investigate «he model's character- 
tecaes ign the sealw of Exman dynmamecs. The curl of she waad 
stress could also be calculated to determine areas of Ekman 
Sice2onmend pumping. Pius, Ge®plex <ransfer functions could 


give an indication of their frequency dapendence. 
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The high frequency peak found in the spectral results 
was unexpected. It is located at a period roughly 2.5 times 
the time step in the fine mesh model and displays vigorous 
eyclonic Wotion. In the coarse grid resolution, it is 
shifted to approximately three times the time step. Its 
amplitude is magnified in the coarser resolution configura- 
ON. A possible explanation is that it is related to the 
"time splitting" problems inharent with the leapfrog time 
differencing scheme. To test this hypothesis, an Euler- 
backward procedure could be inserted every time step to damp 
Out this “time splitting" in the numerical solution. This 
would remove this peak from the spectrum. ie Showed be 
nowed that the energy level associated with this time split- 
ting is very smali in comparison with the currents simulated 
=n the model. 

Further investigat*tion is r2quired to fully understand 
this model's circulation characteristics. Once analysis of 
She N@RPAX d=ifter data 1S complete, a verification of the 
model's performance should be conducted. Thre“ cesuilts Of 
these studies will validate the mcedel and make it a viable 
“ool in understanding the structure and time-varying general 


Gipew@icemon Of the North Pacific. 


65 











APPENDIX A 


The following tables show the results of the rotary 
spectrum analysis performed on the 1976-77 fine mesh data 
set. The model and theoretical values iisplayed in these 
results are from a ratio of the cyclonic to anticyclonic 
components in the motion analyzed. The larger the value 
listed, the stronger the cyclonic motion at that particular 
frequency. Recall from theory that a value of zéro equates 
to the inertial frequency while a value of one corresponds 
to rectilinear motion. the variables pr2sented in the fol- 
lowing tables correspend to: 


FREQ = frequency (cycles/da y) 


Q 
(om 
om 
i) 


Uu-CULrTent autospectrun 

GVV = v-current autospectrum 

QUV = quadrature or inaginary cross spectrum 
of u and v 

MODEL = model rotary spectrum (¢€yuation 2) 


THEORY = thecretical rotary spectrum (equation 3) 
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TABLE IV 


@a 


m), Fine Mesh Da‘ 


$F 43," 


QUV 


Rotary Spectrum k 


Peon Y 


MO DEL 


GVV 


GUU 
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